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Abstract 

A total of 16 BnaGLNI genes coding for cytosolic glutamine synthetase isoforms (EC 6.3.1.2.) were found in the 
Brassica napus genome. The total number of BnaGLNI genes, their phylogenetic relationships, and genetic locations 
are in agreement with the evolutionary history of Brassica species. Two BnaGLNI. 1, two BnaGLNI. 2, six BnaGLNI. 3, 
four BnaGLNI A, and two BnaGLN1.5 genes were found and named according to the standardized nomenclature 
for the Brassica genus. Gene expression showed conserved responses to nitrogen availability and leaf senescence 
among the Brassiceae tribe. The BnaGLNI. 1 and BnaGLNI. 4 families are overexpressed during leaf senescence and 
in response to nitrogen limitation. The BnaGLNI. 2 family is up-regulated under high nitrogen regimes. The members 
of the BnaGLN1.3 family are not affected by nitrogen availability and are more expressed in stems than in leaves. 
Expression of the two BnaGLNI. 5 genes is almost undetectable in vegetative tissues. Regulations arising from plant 
interactions with their environment (such as nitrogen resources), final architecture, and therefore sink-source rela- 
tions in planta, seem to be globally conserved between Arabidopsis and B. napus. Similarities of the coding sequence 
(CDS) and protein sequences, expression profiles, response to nitrogen availability, and ageing suggest that the roles 
of the different GLN1 families have been conserved among the Brassiceae tribe. These findings are encouraging the 
transfer of knowledge from the Arabidopsis model plant to the B. napus crop plant. They are of special interest when 
considering the role of glutamine synthetase in crop yield and grain quality in maize and wheat. 

Keywords: Alloploidization, Brassica napus, Brassica oleracea, Brassica rapa, nitrogen metabolism, nitrogen remobilization, 
senescence. 
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Introduction 

Winter oilseed rape {Brassica napus L.) is the dominant oil- 
seed crop in northern Europe, and nitrogen (N) fertilization 
is the main operational cost for farmers (50% of the total cost 
of production). When compared with other crops, oilseed 
rape is characterized by low nitrogen use efficiency (NUE) 
(Rathke et al, 2006). Despite a high N-uptake efficiency 
(Laine et al. , 1993), only half the N originating from fertilizer 
application is recovered in the seeds (Schjoerring et al, 1995). 
Oilseed rape is characterized by early leaf shedding and unu- 
sual high N loss in senescing falling leaves. The plant can lose 
up to 15% of its entire N content in this way (Rossato et al, 
2001). Leaf senescence generally corresponds to the mobili- 
zation of N reserves from source leaves to sink organs such 
as seeds (Masclaux-Daubresse et al, 2008). In oilseed rape, 
it has been shown that N can be remobilized from senescing 
leaves to expanding leaves at the vegetative stage (sequential 
senescence) as well as from senescing leaves to seeds at the 
reproductive stage (monocarpic senescence) (Malagoli et al , 
2005). 

The rate of senescence and remobilization of leaf N are 
related to the N nutrition status of the plant and to source- 
sink relations (Masclaux et al, 2000). N remobilization 
towards new developing organs is largely dependent on senes- 
cence-related catabolism events and translocation of leaf N 
compounds. Amino acids derived from protein catabolism 
are exported via the phloem to growing parts of the plant; 
the concentration of amino acids in the phloem sap increases 
during leaf senescence (Herrera-Rodriguez et al, 2006, 2007; 
Masclaux-Daubresse et al, 2006). In many species includ- 
ing B. napus, aspartate, glutamate, and their corresponding 
amides are the principal forms of amino N compounds trans- 
ported in the phloem and play a key role in rendering N avail- 
able for remobilization from senescing leaves (Tilsner et al , 
2005). Enzymes involved in the biosynthesis and metabolism 
of amino acids destined for phloem loading are of special 
interest. 

In plants, glutamine synthetase (GS; EC 6.3.1.2) is a key 
enzyme which catalyses an ATP-dependent conversion of 
glutamate to glutamine using ammonium derived from pri- 
mary N uptake and various internal N recycling pathways 
including catabolic release of ammonium during senescence 
(Bernard and Habash, 2009). In a large variety of plants, 
induction of cytosolic glutamine synthetase (GS1) genes 
has been detected during leaf senescence, while chloroplastic 
synthetase isoenzyme (GS2) expression decreases (Masclaux 
et al , 2000; Guo et al , 2004; Martin et al , 2006). It has been 
proposed that in young photosynthetic leaves, the chloroplas- 
tic isoenzyme GS2 is mainly involved in the assimilation of 
ammonium provided by nitrate reduction and photorespira- 
tion through the GS/GOGAT cycle (Masclaux et al, 2001). 
In old senescing leaves, as chloroplasts are breaking down, 
glutamine to be exported would be synthesized by the newly 
expressed cytosolic GS1 isoforms (Masclaux-Daubresse 
et al, 2006). 

The importance of GS1 in N management, growth rate, 
leaf senescence onset and severity, yield, and grain filling 



has been confirmed by co-location of quantitative trait loci 
(QTLs) and functional genomics approaches mainly per- 
formed on maize (Hirel et al, 2001; Martin et al, 2006) and 
rice (Tabuchi et al, 2005). In maize, Glnl.4 is up-regulated 
during senescence (Martin et al, 2005). The Glnl.4 knock- 
out mutation led to a dysfunction in N export and a sharp 
reduction of kernel yield (Martin et al, 2006). GLN1.4 was 
proposed to be involved in re-assimilation of ammonium 
released during leaf protein degradation. In rice, mutants 
lacking OsGSl;l are severely impaired in growth rate and 
grain filling, and glutamine levels in mutant leaf blades are 
reduced (Tabuchi et al, 2005). As the gene product is located 
in companion cells and parenchyma cells of leaf tissues, it has 
been proposed that OsGSl;l is responsible for generation of 
glutamine for remobilization via the phloem. 

To date, all studies on plant genomes have revealed multi- 
genic families coding for several GS1 isoforms. In rice, three 
GLN1 genes have been identified, with seven in wheat, five in 
maize, and five in Arabidopsis thaliana. Transcript omic data 
showed that three A. thaliana genes, AtGLNl.l, AtGLN1.2, 
and AtGLNl.4, are induced during leaf ageing (Guo et al, 
2004). Promoter: :GFP (green fluoresacent protein) fusions 
were used to investigate AtGLNl gene expression in roots. 
AtGLNl.l was localized at the root surface layer, whereas 
AtGLNl. 2, AtGLNl. 3, and AtGLNl. 4 were expressed in root 
vascular tissues (Ishiyama et al, 2004). Detailed expression 
of AtGLNl in leaves was only reported for AtGLNl. 2 that is 
induced in root and leaf tissues under a high N regime and 
is mainly expressed in veins and mesophyll cells in older leaf 
tissues (Lothier et al. , 201 1). In veins, AtGLNl .2 protein was 
localized in the companion cells. The knock-out mutant phe- 
notype led to the conclusion that AtGLNl. 2 is essential for N 
assimilation under ample nitrate supply and for ammonium 
detoxification (Lothier et al, 201 1). For all plant species, it is 
clear that not all GS1 isoforms participate equally in N man- 
agement and remobilization. Regulation of expression is then 
a key clue towards the identification of GLN1 genes poten- 
tially involved in N remobilization. 

Accumulation of GS1 and a decrease in GS2 polypeptides 
were observed in B. napus leaves after onset of leaf senes- 
cence (Ochs et al, 1999). Up to now, four closely related 
genes coding for GS1 isoenzymes, BnGSRl-1, BnGSRl- 
2, BnGSR2-l, and BnGSR2-2, have been identified using 
B. napus root-derived expressed sequence tag (EST) libraries 
(Ochs et al, 1999). Analysis of different tissue types has also 
revealed that these genes are expressed in senescing leaves 
(Buchanan- Wollaston and Ainsworth, 1997). Recent studies 
of Brassicaceae genomes show that the genome of B. napus, 
which is a recent allotetraloid (2^=4x=38, AACC) arising 
from the natural hybridization of monogenomic diploids 
Brassica rap a (A A) and Brassica oleracea (CC) (Nagaharu, 
1935), contains additional genes coding for GS1 isoenzymes. 
Analysis of Brassica lineage genomes revealed that a whole- 
genome triplication occurred shortly after their divergence 
from Arabidopsis (Parkin et al, 2005). Therefore, gene fami- 
lies are more frequent, larger, and more complex in B. napus 
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than in A. thaliana. Brassicaceae genome sequences are also 
highly conserved and many synthenic regions have been iden- 
tified (Paterson et al , 2001; Parkin et al. , 2005; Schranz et al. , 
2006), allowing the identification of 'true' orthologous genes 
between A. thaliana and B. napus. 

In the present study, advantage is taken of the Brassicaceae 
genome structure and of its recent sequencing (unpublished) 
in order to identify all BnaGLNl genes coding for GS1 iso- 
enzymes. It is demonstrated that they are differentially regu- 
lated depending on tissue type, senescence, and N availability. 
The potential role of the BnaGLNl genes in N remobilization 
during leaf senescence, the impact of whole-genome duplica- 
tions and merging on the evolution of the GLN1 multigenic 
family in the Brassiceae tribe, and strategies based on knowl- 
edge transfer from A. thaliana to crop plants are discussed. 

Materials and methods 

Brassica gene identification 

Sequence searches by similarity to A. thaliana AtGLNl coding 
sequences were performed in the GenBank and Genoplante data- 
bases using the BLAST algorithm (Altschul and Lipman, 1990) and 
the A. thaliana AtGLNl coding sequences AtGLNl. 1 (NM_1231 19, 
At5g37600), AtGLNl. 2 (NM_105291, Atlg66200), AtGLNl. 3 
(NM_1 12663, At3gl7820), AtGLNl. 4 (NM_121663, At5gl6570), 
and AtGLNl. 5 (NM_103743, Atlg48470). The BlastN option was 
used to recover B. napus, B. rapa, and B. oleracea ESTs and complete 
mRNAs as well as genomic sequences of B. napus and B. rapa (Cheng 
et al. , 201 1; http://brassicadb.org/brad/). Contigs of ESTs were built 
using the CAP3 algorithm (Huang and Madan, 1999) and validated 
through multiple sequence alignments with ESTs and AtGLNl cod- 
ing sequences, allowing a manual proofreading. Alignments were 
generated with the ClustalW algorithm (Thompson et al, 1994) 
using Clustalw-sequence and Clustalw-profil options available at 
the MOBYLE platform (http://mobyle.pasteur.fr). ESTs included in 
each contig are listed in Supplementary Data File 1 available at JXB 
online. 

BnaGLNl contigs were then enriched and/or their coding 
sequence completed with new cDNA sequences: clones from 
Genoplante oilseed rape cDNA libraries and the ADIZ-MPIZ 
021 library corresponding to ESTs were sequenced when available 
(Supplementary Table SI at JXB online). When the coding sequence 
from cDNAs was incomplete or no clone was available, specific 
primers were designed to clone the total or missing coding region 
(Supplementary Table S2). The amplified fragments were cloned 
into pGEM®-T Easy plasmids (Promega) according to the recom- 
mendation of the supplier, and sequenced. Universal T6 and SP7 
primers, as well as specific primers were used to sequence the clones 
on the positive and negative strands (Supplementary Table S3). All 
DNA sequencing was performed by Cogenics (Grenoble, France) 
and sequences were submitted to GenBank (accession numbers are 
given in Supplementary Table SI). 

Sequence analysis 

A global alignment of coding sequences from mRNA, inferred cod- 
ing sequences, or newly created contigs from ESTs was generated 
with ClustalW (Thompson et al, 1994). Distance matrixes were 
computed using the Dnadist algorithm with a Kimura 2 nucleotide 
substitution model, and bootstrap analysis was performed with 
1000 iterations. A consensus unrooted tree was then generated using 
the Neighbor- Joining method. All algorithms are contained in the 
Phylip 3.67 package available at the MOBYLE platform (http://mob- 
yle.pasteur.fr). The NCBI Conserved Domain Database (Marchler- 
Bauer et al, 2011) was searched with translated B. napus and 



A. thaliana GLN1 -coding sequences. A multiple protein sequence 
alignment was generated with the ClustalW algorithm. 

Genetic mapping and genome or chromosome assignment 

Genetic mapping of BnaGLNl genes was realized using three differ- 
ent B. napus double haploid (DH) populations. The Stellar xDrakkar 
(SD), DarmorxSamourai (DS), and Darmor-Z?z/zxYudal (DY) 
populations consist, respectively, of 94, 134, and 445 genotype DH 
lines described by Lombard and Delourme (2001) and Delourme 
et al. (2006). Gene-specific primers were designed and selected for 
a presence/absence polymorphism in one of the three populations 
(Supplementary Table S2 at JXB online). Linkage analyses were 
performed as previously described by Auger et al. (2009) using 
MAPMAKER/EXP version 3.0b (Lander et al., 1987) and frame- 
work maps constructed in Lombard and Delourme (2001) and 
updated in Delourme et al. (2006). BnaGLNl genes were assigned to 
a linkage group using the ASSIGN command (LOD threshold=8.0) 
and then placed in the most confident interval with the PLACE 
command (LOD threshold=2.0). Recombination frequencies were 
converted into centiMorgans (cM) with the Kosambi function 
(Kosambi, 1944). 

BnaGLNl gene assignment to A or C Brassica genomes was 
performed using a panel of diverse B. napus, B. oleracea, and 

B. rapa genotypes available in the authors' group. The panel con- 
tains genomic DNA from B. napus genotypes Darmor-/?z/z, Yudal, 
Stellar, Drakkar, Samourai, aburamassari, Aviso, Tenor, Express, 
Montego; B. rapa genotypes Zl, CI. 3, Chiifu; and B. oleracea geno- 
types HDEM and CI 02. This panel was PCR screened with specific 
but not polymorphic gene markers (Supplementary Table S2 at JXB 
online). 

Chromosome assignment for BnaC.GLNl genes was realized 
using monosomic and polysomic addition lines carrying one or 
several additional C chromosomes from Darmor-/?z/z. Lines were 
selected from a cross between B. napus Darmor-Z?z/z and B. olera- 
cea CI. 3 (A.M. Chevre and F. Eber, INRA Rennes, unpublished 
results). Genomic DNA from these lines was PCR screened with 
specific but not polymorphic markers (Supplementary Table S2 at 
JXB online). 

Plant material and growth conditions 

Brassica napus L. plants from the Darmor-/?z/z genotype were grown 
in a greenhouse at INRA Versailles, France. Seeds were sown on 
sand and watered with nutritive solution during 2 weeks in order to 
allow germination and subsequent growth of plantlets. When the 
first two true leaves appeared, plantlets were transferred into pots 
containing sand and were separated into two groups with contrast- 
ing N fertilization regimes (LN for low nitrate and HN for high 
ntrate, 0.4 mM and 8mM N0 3 ", respectively) according to Albert 
et al. (2012). At 56 d after sowing, four plants of each nutrition 
regime were harvested and sampled. For each plant, all leaf ranks 
were collected: primary and secondary veins were separated from 
the rest of the leaf, described as the limb. All fresh samples were 
frozen immediately in liquid nitrogen and stored at -80 °C. 

Brassica napus L. plants from the Express genotype were grown 
in field trials in 2009-2010, in Le Rheu (Brittany), France. Seeds 
were sown on 7 September 2009 with plant density set at 40 plants 
m~ 2 . The field trial was conducted with contrasting N fertilization 
regimes. Plant N status was monitored over the vegetative stage 
by calculating the nitrogen nutrition index (NNI) (Colnenne et al. , 
1998). The balance-sheet method was used as a decision tool for N 
fertilization, setting the potential yield for LN and HN regimes at 20 
q ha" 1 and 35 q ha" 1 , respectively (Makowski et al. , 2005). LN plants 
did not receive N fertilizer, while the HN plants received a total input 
of 110 kg N ha" 1 spread at two different times (12 February and 19 
March 2010). LN and HN plants were harvested, respectively, on 9 
and 12 April, at the beginning of the flowering period when half the 
plants of the plot had their first flowers open on the main stem (Fl, 
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or 60 on the BBCH scale), and 400 degree-days later (base 0) on 17 
and 20 May at the beginning of the seed filling period (G2, or 71-73 
on the BBCH scale). 

Plants from 0.5 m 2 per plot (~20 plants) were harvested in the 
early morning and sampled during the subsequent hour. For each 
batch, plants were ranked according to their length and develop- 
mental stage; the six median plants were selected for sampling. On 
the main stem, the lowest leaf starting to yellow (Old) and the high- 
est leaf at least 5 cm long (Young) were selected, the petiole and 
main vein were removed, and limbs were sampled. The stems above 
young and old leaf insertions were selected and sampled over 2 cm 
and 4 cm, respectively. All fresh samples were frozen immediately in 
liquid nitrogen and stored at -80 °C. 

Nucleic acid manipulation 

PCRs were conducted in a 20 ul mix containing 2-10ng of DNA, 
0.25 mM dNTPs (Promega), 0.5 uM of each primer (Eurogentec, 
Angers, France), and 0.5 U of Taq DNA polymerase (Promega) 
in the appropriate buffer supplemented with 2.5 mM MgCl 2 . The 
amplification program was run on a PTC-225 thermocycler (MJ 
Research, Waltham, MA, USA) with the following conditions: 35 
cycles of denaturation at 94 °C for 30 s (3min for the first cycle), 
annealing at 55-60 °C for 30 s, and elongation at 72 °C for 1-2 min 
(lOmin for the last cycle). 

Total RNAs were extracted with the SV Total RNA Isolation 
System (Promega) from 70 mg fresh weight (FW) of ground frozen 
tissue. First, samples were homogenized in the RNA lysis buffer 
(400 ul) using TissueLyserll from Qiagen. Then, all cell debris was 
eliminated by filtering the lysate through a 'Nucleospin 96 RNA fil- 
ter Plate' (Macherey-Nagel). The manufacturer's protocol was then 
followed. In order to remove any remaining DNA traces, 1.5 jug 
of RNA was treated with DNase using the Turbo DNA-free kit 
(Ambion) according to the manufacturer's protocol. The quality 
of RNA was assessed by an electrophoresis on agarose gel (1.3%, 
w/v), and the absence of DNA contamination in samples was con- 
firmed by PCR amplification. First-strand cDNA was synthesized 
using 2 |ug of total RNA with oligo(dT) 12 _i 8 primers and Superscript 
III reverse transcriptase (Invitrogen) according to the manufactur- 
er's instructions. Two independent reverse transcription reactions 
were performed as technical replicates. cDNA samples were diluted 
26-fold with sterile water before use. 



Gene expression analysis with qPCR 

Quantitative PCRs (qPCRs) were set up with the LightCycler 480 
SYBR Green I Master mix (Roche Diagnostics) and 4 ul of diluted 
cDNA in a final volume of 12 ul. The concentration of specific for- 
ward and reverse primers was set at 0.42 uM. qPCRs were run on a 
Light Cycler LC480 (Roche Diagnostics) under the following condi- 
tions: an initial step at 95 °C for 10 min, then 50 cycles of 95 °C for 
10 s, 60 °C for 10 s, and 72 °C for 20 s. 

Primer pairs were designed for short and specific amplification of 
individual members of the multigenic BnaGLNl family and five ref- 
erence genes (Supplementary Table S4 at JXB online). Identification 
of B. napus reference genes for reverse transcription-qPCR analy- 
sis was based on EST sequence similarity to A. thaliana genes with 
verified expression stability over a wide range of tissues and growing 
conditions. Arabidopsis thaliana genes were selected from among the 
list established by Czechowski et al. (2005). Their respective coding 
sequences were used to retrieve highly similar B. napus ESTs from 
GenBank using the BlastN algorithm. One B. napus EST per ref- 
erence gene was selected to design qPCR primers (Supplementary 
Table S4). 

For each run, single product amplification was confirmed by melt 
curve analysis. PCR products from each primer pair and genotype 
were sequenced in a preliminary analysis. The amplification effi- 
ciency was assessed for each genotype, with each primer pair using 
a dilution curve method over six orders of magnitude, on a pool of 



cDNAs from different tissues and modalities. Selected primer pairs 
have efficiencies >1.8. 

The results reported were obtained from four biological replicates 
and two reverse transcriptions as technical replicates. All samples, 
including reverse transcription and biological replicates, were run 
at the same time for each primer pair. Raw fluorescence data were 
collected and analysed with the R package 'qpcR' (Ritz and Spiess, 
2008). The 'pcrbatch' function was used to select sigmoid models 
for the fluorescence curves and then allowing the determination of 
the intrinsic amplification efficiency (sig.eff) and threshold cycle 
(sig.CpD2) at the second derivative maximum (Rutledge, 2004). For 
each run, cDNA relative quantity (RQ) was calculated using the effi- 
ciency mean value from the two technical and the four biological 
repetitions (mean sig.eff), and the run-specific threshold cycle as: 
RQ=l/(mean sig.eff) (sig.CpD2). 

The most stable reference genes were selected using the GeNorm 
method from Vandesompele et al. (2002) available through the 
'SLqPCR' R package (Dr Matthias Kohl SIRS-Lab GmbH). The 
four reference genes BnaX.PTB, BnaX. SAND, BnaXPP2A, and 
BnaX. UBC21 from the five tested were retained with an average M 
value equal to 0.33. For each cDNA sample, a normalization factor 
(NF) was calculated as the geometrical mean of RQ from the four 
selected genes, and normalized RQ (NRQ) was then calculated as 
NRQ=RQ/NF. The mean of both technical replicates was then cal- 
culated for each sample. 



Results 

Identification ofGLNI coding sequences in EST 
databases of Brassica napus and its progenitors 
Brassica oleracea and Brassica rapa 

A total of 588 B. napus ESTs, 126 B. rapa ESTs, and 36 
B. oleracea ESTs highly similar to one or several of the 
five AtGLNl.l-AtGLN1.5 mRNA-coding sequences from 

A. thaliana were isolated from public and private Genoplante 
databases (Supplementary Data File SI at JXB online). 
Sequence assembly and alignments with AtGLNl mRNA 
sequences revealed different groups of transcripts that 
allowed ESTs to be grouped and 16 individual contigs that 
might correspond to different BnaGLNl genes of B. napus 
to be extracted (Supplementary Data File S2). Eight con- 
tigs for B. rapa and seven contigs for B. oleracea were also 
isolated (Supplementary Data Files SI, S2). The BnaGLNl, 
BraGLNl, and BolGLNl contigs from B. napus, B. rapa, and 

B. oleracea, respectively, show high levels of sequence simi- 
larity with all the five AtGLNl genes (Table 1). The analysis 
of similarity levels and a phylogenetic tree (Fig. 1) allowed 
homologous sequences for each AtGLNl mRNA to be 
clearly identified in B. napus, B. rapa, and B. oleracea. The 
phylogenetic tree reveals that the GLN1 mRNA sequences 
are divided into two distinct groups for monocotyledonous 
(wheat, maize, and rice) and dicotyledonous species. The 
Brassicaceae sequences are divided into five clusters, each 
one including one A. thaliana and one or several B. napus, 
B. oleracea, and B. rapa sequences (Fig. 1). Each B. napus 
sequence is closely related to one sequence from either 
B. oleracea or B. rapa progenitors, illustrating the ances- 
tral relationship with the A and C Brassica genomes. Each 
B. napus sequence is also closely related to one of the five 
A. thaliana AtGLNl mRNA sequences, allowing the iden- 
tification of homeologous related sequences between the 
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Sequence name Gene name 



-HZ 




8. rapa 
8. oleracea 
B. napus 
A. thaliana 



nDUOf \J\3\J 


UouO / / 




ukm i ztzuy 


/ duo / d 




UW I I U 


Tad^lh 
/ duo / JU 




nP)194911 


/ duO / O 




VfiKQOQ 


7mdl Kid - d^1 

Z-lll\Jl-l\l*t — DO / 




YfiCQOQ 

Aooyzy 


/.1/luL/vO — DO / 




nD I OUDOC7 


UouO / O 




Xfi^997 


7mdl NI9 - d^i 

Z.///Oi_/VZ — uo / 


2 


AY49197D 


/ duOc / 




AY4Q1Q71 
nl ty I c? f I 


ToOCoO 
/ duoct 




aalj I ouuoo 






AY4919F58 


TaGSrl 




AY4Q1QfiQ 


TaGSr2 




Xfi^99fi 


7mdl A/fi = d^1 

Z.///UL/VU 00 / 


•/ 


INIVI UU I I I lOZf 


7mdl A/9 - f^Q-/ 
Z./;/OL/Vv5 — 00 / 


5 


NM 10^74^ 
in ivi i y/tj i *ttj 


f~\ IWL/V 1 . 




RnlPlnl *=> P1 a 

DUIOlll I .D u I 


RnIC dl KM ^ a 




DlldOIII 1 .O u 1 


RnaC HI KM ^ a 
Di WO. OL/V / .O.d 




RnaGlnl ^ P9 


RnaA CSI KI1 ^ a 




RraGlnl-'SPI 3 

□ ICIVJMI 1 Ju 1 


RmA CSI KI1 f> 

Dl dA\. uLl v / .O.d 


< 

* 


INIVI 1 l£UUJ 


AtGI N1 ? 

n(wL/V f .O 


RnlPlnl ^ P9 a 
DUIOlll 1 .O 


DU/u. Oi_/V / . v5. 0 




RraPlnl ^ P9a 
DldOIII 1 .O 


Rra4 ^1/ A/-/ 9 r 




RnaGlnl ^ PR 

Dl IdOII 1 1 .O \j\J 


RnaC HI KM 9 r 




J/\OUUUC7 i + 


RnaA dl KM 9 r 




RnaGlnl ^ P^ 

Dl IdVJII 1 1 .O UO 


RnaC Cll KI1 9 h 




RnaGlnl ^ P4 

Dl IdOII 1 1 .O 


RnaA HI KM 9 h 




RraGlnl 1 P^ 

Dl avJll 1 1 . 0 UO 


RraA dl KI1 In 




JAOUDDSO 


RnaA HI KM 9 a 




RraGlnl ? P1 

Dl dOII I I . O u I 


RraA dl KI1 9 a 

D/arl. Oi_/V 1 .O.d 




O /\OU\J\J<j\j 


RnaC dl KI1 9 a 
Dllau.uLlV I .O.d 




RnlPlnl ^ P1 
DUIOlll I .O u I 


Rnld dl KM 9 £) 

DU/u. Oi_/V 1 .O.d 




NM 191fifi^ 

INIVI 1 £ 1 UUO 


AtGLN1.4 




IY?nfifiQ7 

JAOUOD3 1 


RnaA dl KM d a 




RraGlnl 4 P1 

Dl dOII I I .t U I 


RmA dl KI1 d a 




U /\ \J \J \J \J \J \J 


RnaC, dl N1 d a 




RnlGlnl 4 P1 

DU I \J 1 1 1 1 .t U | 


Rn/P A/"/ 4 p 




BnaGlnl 4 C4 


Rna4 dl M1 d h 

LDI /C?/i . UL / V / .'-r.kJ 




RraGlnl 4 P9 

Dl dOII 1 1 .t 


RraA dl KM d h 

Dl dry. OZ_/V / .H.U 




RnaGlnl 4 C\ 

Dl IdOII 1 1 wO 


RnaC dl KI1 d h 

Dl IdO. OL/V / .H.U 




BolGlnl 4 G2 


RnIC dl N1 d h 




NM 19^11Q 

IN IVI 1 1 1 9 


AtGLN1.1 


X82997 CBnGSR2-'h 


RnaA dl M1 1 a 




RraGlnl 1 P1 

Dl dull I I . I u I 


Rra4 dl KI1 1a 
di dr\. yji—i v / . / . d 




BnaGlnl. 1 C2 


BnaC.GLN1.1.a 




DUioin i . i_u i 


dO/u.ul/v i . i .a 


-< 


NM_1 05291 


AtGLN1.2 


X76736 (BnGSRM) 


BnaA.GLN1.2.a 




AY773089 


BraA.GLN1.2.a 




Y12459 (BnGSR1-2) 


BnaC. GLN1.2.a 




EU822334 


BolC.GLN1.2.a 





Fig. 1. Cytosolic glutamine synthetase (GS1) phylogenetic tree. DNA coding sequences (CDS) were aligned using Clustal. The distance matrix 
was computed using Dnadist with a Kimura 2 nucleotide substitution model (bootstrap analysis, 1000 iterations). A consensus unrooted tree was 
then generated using the Neighbor-Joining method from the Phylip 3.67 package. Black and grey dots indicate bootstrap values >90% and 50%, 
respectively. All programs are available at www.mobyle.pasteur.fr. incomplete CDS sequence when compared with the A thaliana reference sequence. 



four species. The contigs were then named according to the 
AtGLNl.x (x from 1-5) gene with the highest sequence simi- 
larity. Since for each AtGLNl.x sequence several BnaGLNl 
sequences were found, contig names were also extended by 
a copy number (Cn): BnaGLNl. x_Cn (Tables 1, 12; Fig. 1). 
Similarly, the names of the B. rapa and B. oleracea contigs 
follow the same rules (BraGLNl.x_Cn and BolGLNl.x_Cn, 
respectively; Table 2). 



Most of the EST assemblies have been confirmed through 
sequencing partial or full-length cDNA clones when available 
(Supplementary Table S 1 at JXB online). For the few EST assem- 
blies that cannot be confirmed in this way, the cloning of missing 
coding sequences (CDS) was performed by designing primers 
from the B. rapa and B. oleracea homologous contig sequences 
(see the Materials and methods). This allowed the completion of 
the BnaGLNl. 4_C3 and BnaGLNl. 4_C4 sequences. 
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Table 2. Names ofcontigs ofESTs, mRNA, and AtGLNI homologous genes in Brassica napus, B. oleracea, and B. rapa 

The names of genes encoding each contig were assigned according to the Ostergaard and King (2008) nomenclature, taking into account A or 
C genome location, and the closest AtGLNI sequence homology. 



A. thaliana 


Brassica napus 




Brassica oleracea 




Brassica rapa 




Gene name 


Name of contig of 


Gene name 


Name of contig of 


Gene name 


Name of contig of 


Gene name 




ESTs; mRNA name 




ESTs; mRNA name 




ESTs; mRNA name 


(BRAD name; LG) a 


AtGLNI. 1 


BnaGLN1.1_C1;X82997 


BnaA.GLNI . 1.a 






BraGLN1.1_C1 


BraA.GLNI .1 .a 


(At5g37600) 


(BnGSR2. 1) 










(Bra028132; A04) 




bnaULNi .\_L>z, Yiz4bU 


Dnau.LiLN I . i .a 


Dz-nI^i km -\ r^-\ 
bOlULNl 


dOIL>.kjLN 1 . 1 .a 








(BnGSR2.2)° 












AtGLNI. 2 


BnaGLN1.2_C1;X76736 


BnaA.GLNI. 2.a 






BraGLN1.2_C1; 


BraA.GLNI. 2.a 


(At1g66200) 


(BnGSR1.1) 








EU499383; AY773089 


(Bra039756; A02) 




BnaGLN1.2_C2; Y12459 


BnaC.GLN1.2.a 


BolGLN1.2_C1; 


BolC.GLN1.2.a 








(BnGSRI .2) 




EU822334; EU822335 








AtGI N1 9 


RnaGI N1 9 OP- IX90RRQ9 


RnaA Gl N1 In 






RraGI N1 S 01 

LJI CL\J L_l N I .O w I 


RraA Gl N1 S a 

LJ I d/A . \J l_l N I . O . CI 


(At3g 17820) 












(Bra022247; A05) 




BnaGLNI .3_C1 ; JX306690 


BnaC.GLN1.3.a 


BolGLN1.3_C1 


BolC.GLN1.3.a 








BnaGLN1.3_C4 


BnaA.GLNI. 3.b 






BraGLN1.3_C3 


BraA.GLNI. 3.b 

(Bra021276; A01) 




BnaGLNI. 3_C3 


BnaC.GLN1.3.b 












BnaGLNI .3_C5; JX306694 


BnaA.GLNI. 3.c 






BraGLNI .3_C2 b 


BraA.GLNI. 3.c; 
(Bra001686; A03) c 




BnaGLNI .3_G6 


BnaCGLNi .3.c 


BolGLNI .3_G2 


B0ICGLN1 .3.C 






AtGLNI. 4 


BnaGLNI .4_C1 ; JX306697; 


BnaA.GLNI A.a 






BraGLN1.4_C1 


BraA.GLNI. 4.a 


(At5g 16570) 


JX306692 C 










(Bra023573; A02) 




BnaGLNI .4_C2; JX306695 


BnaC.GLN1.4.a 


BolGLN1.4_C1 


BolC.GLN1.4.a 








BnaGLNI .4_C4; JX306700 C 


BnaA.GLNI A.b 






BraGLNI. 4_C2 


BraA.GLNI. 4. b 

(Bra008612; A10) 




BnaGLNI .4_C3; JX306698 C 


BnaC.GLN1.4.b 


BolGLN1.4_C2 


BolC.GLN1.4.b 






AtGLNI. 5 


BnaGLNI. 5_C2 


BnaA.GLNI. 5.a 






BraGLNI. 5_C1 b 


BraA.GLNI. 5.a 


(At1g48470) 












(Bra018729; A06) 




BnaGLNI. 5_C1 ; JX306691 b 


BnaC.GLN! .5.a 


BolGLN1.5_C1 b 


BolC.GLN1.5.a 







a Annotation and localization on the linkage group from BRAD, the Brassica rapa genome sequencing project consortium (Wang etal., 201 1). 
b Incomplete CDS sequence when compared with the A. thaliana CDS. 

c SNP insertion disrupting the ORF when compared with the A. thaliana CDS reference sequence and the BnaGLNI contig. 



Except for the four BnaGLNI. 1_C1, BnaGLNI. 1_C2, 
BnaGLNI. 2_C1, and BnaGLNI. 2_C2 contigs, no other 
BnaGLNI sequences have ever been described previously in 
the literature or reported in databases as glutamine synthetase 
gene products. While the BnaGLNI. 1_C1, BnaGLNI. 1_ 
C2, BnaGLN1.2_Cl, and BnaGLN1.2_C2, sequences have 
been identified as BnGSR2-l, BnGSR2-2, BnGSRl-1, and 
BnGSRl-2 mRNA, respectively (Table 2; Ochs et aL, 1999), con- 
tig analysis allowed the completion of the 5' end (untranslated 
region and coding sequence) of the BnGSR2-2 sequence that was 
previously missing (Supplementary Data File S3 at JXB online). 

Genetic localization of the BnaGLNI loci on the A or 
C Brassica genome using PCR and gene name 
annotation 

Phylogenetic analyses showed a strong relationship between 
each BnaGLNI gene and a gene from one or other of the pro- 
genitors B. rapa and B. oleracea, suggesting a common ancestral 



origin on the A or C Brassica genome, respectively. The phylo- 
genetic tree also shows that each B. napus sequence, related to 
one sequence from either progenitor, is also related to another 
B. napus sequence, itself related to the other progenitor. The 
two B. napus homeologous genes, the B. rapa and the B. oler- 
acea genes, are thus defining in this way a homeology group 
(a, b, or c). It was found therefore that each of the AtGLNI. 1, 
AtGLNI. 2, and AtGLNI. 5 genes is related to one homeology 
group, while the AtGLNI. 3 and AtGLNI. 4 genes are related to 
three and two groups, respectively. It has to be noted that the 
b group related to AtGLNI. 3 is incomplete as no BolGLNI. 3 
expressed sequence has been identified. Both the homeology 
groups and the Brassica genome were used to ascribe names 
to the BnaGLNI genes; thus, the genes are named BnafA or 
C genome] GLN 1.x [a, b, or c homeology group] according to 
Ostergaard and King (2008). A similar notation was used for 
the BraGLNI and BolGLNI genes (Table 2). 

In order to identify the A or C genome origin, the poten- 
tial localization of the BnaGLNI genes on linkage groups 
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and/or chromosomes known to arise from the A or C 
Brassica genomes was then investigated. Specific primer 
pairs were designed to localize each BnaGLNl gene 
(Supplementary Table SI at JXB online). Five genes were 
mapped in this way on at least one of the three mapping 
populations available (Stellar xDrakar, DarmorxSamourai, 
and Darmor-fc/zxYudal), recording the presence/absence 
of polymorphism. In good agreement with the phyloge- 
netic tree analysis, the three genes BnaA.GLNI. 3. c, BnaA. 
GLNU.a, and BnaA.GLN1.5.a were localized on linkage 
groups associated with the A genome on chromosomes A03, 
A05, and A06, respectively. The two genes BnaC.GLN1.2.a 
and BnaC.GLN1.3.c were localized on linkage groups asso- 
ciated with the C genome on chromosomes C02 and C03 
(Table 3; Supplementary Fig. SI at JXB online). 

For the other members of the BnaGLNl gene family, an 
attempt to assign the BnaGLNl genes to the A or C genomes 
using mapping populations was unsuccessful. Therefore, a 
panel of various Brassicaceae genotypes (Supplementary Fig. 
S2 at JXB online) was used in order to detect the BolGLNl 
and BraGLNl orthologous genes, using specific BnaGLNl 
primers (Supplementary Table SI). The number of genotypes 
used for each Brassica species was adjusted in order to take 
into account the possible allelic variations and to detect the 
presence/absence of polymorphism. Furthermore, additional 
lines carrying the full A genome and one or several B. napus 
C chromosomes were used in order to determine preciselg 
the localization of the BnaC. GLN1 genes (Auger et al. , 2009) 
(Supplementary Fig. S3). 

Table 3. Genetic mapping of BnaGLNl genes 



The results are summarized in Table 3 and detailed 
in Supplementary Figs SI, S2, and S3 at JXB online. In 
brief, with the exception of BnaC.GLNI. 4. b, all the identi- 
fied BnaGLNl genes were assigned to the A or C Brassica 
genome, confirming the two by two relationship of homology 
between them, which allowed them to be named according to 
their ancestral genome origin as recommended by Ostergaard 
and King (2008). 

Identification of GLN1 genes of Brassica rapa and 
Brassica napus genomes 

The recently sequenced and annotated B. rapa genome 
(BRAD; Cheng et al., 2011) was used to perform BLAST 
searches and sequence alignments using the BraGLNl contig 
sequences identified here. Analysis revealed eight annotated 
BraGLNl genes (Table 2). Alignments between BraGLNl 
gene sequences and contigs revealed potential splicing vari- 
ants. Indeed, the CDS deduced from the BraGLNl. 3 _C2 con- 
tig appeared incomplete at the 5' end. The most highly similar 
Bra001686 annotated gene on the A03 chromosome also 
appeared incomplete when compared with the AtGLNl.3 
CDS, as it is missing the first exon. The BLAST search on 
A03 chromosome vl.l revealed the presence of a sequence 
highly similar to the AtGLNl.3 first exon, 4kb upstream 
of the identified Bra001686 sequence (bp 17 854 801 to 17 
854 849). According to the BRAD annotation, this inserted 
region has been described as an long terminal repeat (LTR) 
transposon of 3746 bp on the minus strand. The identified 



Mapping was performed using specific primers for each contig and different mapping populations or genotypes for B. napus, B. oleracea, and 
B. rapa. Linkage groups used to assign each gene to A or C genomes are presented. 



BnaGLNl Mapping Linkage group LOD Upstream marker Downstream marker 

population (previous name) 



Gene name 


Contig name 








Name 


Distance (cM) 


Name 


Distance (cM) 


BnaA.GLNI . 1.a 


BnaGLNl. 1_C1 


Panel 3 


A 












BnaC.GLNI . 1.a 


BnaGLNl. 1_C2 


Add b 


C06 












BnaA.GLNI .2. a 


BnaGLNl. 2_C1 


Panel 3 


A 












BnaC.GLNI. 2.a 


BnaGLNl. 2_C2 


SD C 


C02 (SD02) 


13.5 


PFM504 


9.1 


J15.1200 


9.1 


BnaA.GLNI .3.a 


BnaGLNl. 3_C2 


DS d 


A05 (DS19) 


17 


E1M4.21 


2.7 


BN614 


13.5 


BnaC.GLNI .3.a 


BnaGLNl. 3_C1 


Add b 


C05 












BnaA.GLNI .3. b 


BnaGLNl. 3_C4 


Panel 3 


A 












BnaC.GLNI. 3.b 


BnaGLNl. 3_C3 


Add b 


C01 












BnaA.GLNI .3.c 


BnaGLNl. 3_C5 


DS d 


A03 (DS04) 


16.1 


PFM 193 


3.8 


BN466 


7.9 


BnaC.GLNI .3.c 


BnaGLNl. 3_C6 


SD C , add b 


C03 (SD717) 


25.9 


BN04C 


1.3 


IGF0193C 


2.2 


BnaA.GLNI A.a 


BnaGLNl. 4_C1 


Panel 3 


A 












BnaC.GLNI A.a 


BnaGLNl. 4_C2 


Panel 3 


C 












BnaA.GLNI A.b 


BnaGLNl. 4_C4 


Panel 3 


A 












BnaC.GLNI A.b 


BnaGLNl. 4_C3 
















BnaA.GLNI .5. a 


BnaGLNl. 5_C2 


DY e 


A06 (DY06) 


8.9 


nr 


nr 


BN57463 


21.8 


BnaC.GLNI .5.a 


BnaGLNl. 5_C1 


Panel 3 


C 













a Panel of B. napus, B. oleracea, and B. rapa genotypes. 

b Monosomic and polysomic addition lines obtained from the Darmor-jbz/7xC1 .3 cross. 
c StellarxDrakar mapping population. 
d DarmorxSmoura! mapping population. 
e Darmor-jbz/ixYudal mapping population. 



Exploring glutamine synthetase gene expression in Brassica napus L | 3935 



B. rapa EST (EX089134) that allowed identification of the 
5 7 region of the BraGLN1.3_C2 contig starts in the trans- 
poson region and continues into the first identified exon of 
Bra001686 which corresponds to the second AtGLNl.3 exon. 

The sequencing programme performed at URGV allowed 
identification of BnaGLNl sequences in the B. napus Darmor- 
bzh genome (SEQ-POLYNAP, ANR-09-GENM-021). The 
BraGLNl protein sequences were deduced from the iden- 
tified BraGLNl genes using the BRAD tool, and used to 
search the database of protein sequences built from B. napus 
genomic sequence analysis (unpublished data). From the 
BnaGLNl protein sequences identified, genomic sequences 
were recovered (Supplementary Data File S4 at JXB online). 
Interestingly, 16 putative BnaGLNl genes and two putative 
BnaGSL genes (coding for the GS2 isoform) were found. 
The 16 BnaGLNl genomic sequences (Supplementary Data 
File S4) were used to analyse similarities with the B. napus 
contigs and to create a phylogenetic tree (Supplementary 
Fig. S4). Interestingly, each genomic sequence was closely 
associated with one contig sequence, suggesting that all 
the genes with the 16 BnaGLNl contigs had been found. 
The deduced mRNA sequences (Supplementary Data 
File S5), obtained using FGENESH software available on 
the SoftBerry website (http://linuxl.softberry.com/berry. 
phtml?topic=fgenesh&group =programs&subgroup=gfind), 
showed very a high similarity with the contig sequences 
(Table 4) and allowed the gene structures to be deduced 
(Fig. 2). Except for BnaGLNl. 3_C5 and BnaGLNl. 3_C6, 



similarities between mRNA and associated contigs were near 
100% (Table 4). Knowing that contig sequences and mRNA 
sequences are obtained from different B. napus genotypes 
(Supplementary Data File SI), this indicates that there is 
almost no polymorphism between the different BnaGLNl 
coding sequences regarding the various genotypes of B. napus 
used for genome and EST sequencing. The BnaGLNl genes 
contained between seven and 12 exons. GLN1.4 and GLN1.5 
genes have the same number of exons in both B. napus and 
Arabidopsis (Table 4). For the other GLN1 genes, exon num- 
bers are different between Arabidopsis and B. napus, but quite 
close; for example, AtGLN1.3 and BnaGLNl. 3 contain fewer 
exons than other AtGLNl and BnaGLNl genes. 

BnaGLNl protein sequence conservation 

Protein sequences of the BnaGLNl family deduced from 
the coding sequences of contigs or from the deduced mRNA 
are similar. The BnaGLNl proteins share between 93% and 
96.6% identity with the AtGLNl proteins encoded by their 
respective orthologous genes (Table 1). Within each homeol- 
ogy group, the A and C BnaGLNl proteins share 98.3-100% 
identity. 

In all BnaGLNl protein sequences, two conserved pfam 
domains specific to glutamine synthetase enzymes (pfam 
03951 and pfam 00120) were identified (Fig. 3). The resi- 
dues involved in the ammonium/glutamate-binding pocket 
(Eisenberg et al , 2000) are also strictly conserved. In contrast, 



Table 4. Comparison between contigs ofESTs and the mRNA sequences deduced from the BnaGLNl genomic sequences 

Comparison of nucleotide length (bp), % similarities, and exon numbers between the B. napus contigs and the mRNA sequences 
(Supplementary Data File S5 at JXB online) deduced from the BnaGLNl genomic sequences (Supplementary Data File S4). Similarities were 
estimated using BLAST (NCBI) and exons using FGENESH software at the SoftBerry website. 



Gene name Contig name Contig Deduced Putative mRNA % similarity No.of exons 

length (bp) mRNA name length (bp) between contig 

and mRNA 



AtGLNl. 1 








1494 




9 


BnaA.GLNI . 1.a 


BnaGLNl. 1_C1 


1374 


mRNA.BnaA.GLN1 .1 .a 


1746 


100 


11 


BnaC.GLNI . 1.a 


BnaGLNl. 1_C2 


1367 


mRNA.BnaC.GLN1 .1 .a 


1474 


100 


11 


AtGLNl. 2 








1499 




10 


BnaA.GLNI .2. a 


BnaGLNl. 2_C1 


1430 


mRNA.BnaA.GLN1 .2. a 


2587 


100 


9 


BnaC.GLNI. 2. a 


BnaGLNl. 2_C2 


1431 


mRNA.BnaC.GLN1 .2. a 


2123 


99 


11 


AtGLNl. 3 








1341 




9 


BnaA.GLNI .3.a 


BnaGLNl. 3_C2 


1336 


mRNA.BnaA.GLN1 .3. a 


2012 


100 


8 


BnaC.GLNI. 3. a 


BnaGLNl. 3_C1 


1487 


mRNA.BnaC.GLN1 .3. a 


1555 


99 


8 


BnaA.GLNI .3.b 


BnaGLNl. 3_C4 


1273 


mRNA.BnaA.GLN1 .3.b 


2051 


99 


7 


BnaC.GLNI .3.b 


BnaGLNl. 3_C3 


1275 


mRNA.BnaC.GLN1 .3.b 


1803 


100 


7 


BnaA.GLNI .3.c 


BnaGLNl. 3_C5 


1253 


mRNA.BnaA.GLN1. 3. c 


1840 


96 


9 


BnaC.GLNI. 3.c 


BnaGLNl. 3_C6 


1245 


mRNA.BnaC.GLN1 .3.c 


1761 


96 


9 


AtGLNl. 4 








1269 




12 


BnaA.GLNI A.a 


BnaGLNl. 4_C1 


1273 


mRNA.BnaA.GLN1. 4. a 


1461 


100 


12 


BnaC.GLNI A.a 


BnaGLNl. 4_C2 


1259 


mRNA.BnaC.GLN1 .4. a 


1415 


99 


12 


BnaA.GLNI A.b 


BnaGLNl. 4_C4 


1102 


mRNA.BnaA.GLN1. 4. b 


2582 


99 


12 


BnaC.GLNI A. b 


BnaGLNl. 4_C3 


1123 


mRNA.BnaC.GLN1. 4. b 


1494 


100 


12 


AtGLNl. 5 








1307 




10 


BnaA.GLNI .5. a 


BnaGLNl. 5_C2 


1392 


mRNA.BnaA.GLN1. 5. a 


1329 


99 


10 


BnaC.GLNI. 5. a 


BnaGLNl. 5_C1 


1380 


mRNA.BnaC.GLN1 .5. a 


1349 


100 


10 
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Fig. 2. Structure o\BnaGLN1 genes. For each BnaGLNI gene, the length of the 5' and 3' untranslated regions (UTRs) (white boxes), exons (black 
boxes), and introns (black lines) is represented by a number corresponding to base pairs. 
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Fig. 3. Alignment of Brassica napus and Arabidopsis thaliana GS1 proteins. Protein sequences were deduced from DNA coding sequences and aligned 
using Clustal. Stars indicate residues involved in the ammonium/glutamate-binding pocket (Eisenberg etal., 2000). Boxes indicate residues involved in 
ammonium affinity properties (Ishiyama etal., 2006). Arrows indicate conserved domains (1) pfam 03951 Gln-synt_N glutamine synthetase bet-Gasp 
domain; and (2) pfam 00120 gln-synt_C catalytic domain. Residues are coloured according to their polarity properties (neutral non-polar as black, neutral 
polar as green, acidic as red, and basic as blue). 



the polar amino acids Q49 and SI 74, shown to be involved 
in the ammonium high affinity properties of AtGLNl.l and 
AtGLN1.4 (Ishiyama et al, 2006), are not strictly conserved 
in all the BnaGLNl.l and BnaGLN1.4 proteins. The polar 
Q49 was converted into an acidic glutamate E49 in all the 
BnaGLNl . 1 and BnaGLNl .4 sequences, and the SI 74 is con- 
served only in the two BnaGLNl. 4. b sequences but was con- 
verted into an A174 in the BnaGLNl. 4.a and BnaGLNl.l. a 
sequences. Depending on the effect of such amino acid modi- 
fications, it might be possible that ammonium affinity prop- 
erties have not been conserved within the BnaGLNl.l and 
BnaGLNl. 4 protein families. In contrast, the residues K49 
and A 174 present in the low affinity enzymes AtGLN1.2 
and AtGLN1.3 are conserved in all the BnaGLNl. 2 and 
BnaGLNl. 3 protein sequences, suggesting the conservation 
of the low ammonium affinity properties in those two protein 
families (Fig. 3). 

Expression of BnaGLNl genes is modified depending 
on the nitrogen regime and leaf senescence 

A first analysis of EST distribution between libraries and 
BnaGLNl contigs led to the conclusion that BnaGLNl genes 



are probably differentially expressed according to tissue and 
developmental stage (Supplementary Data File SI at JXB 
online). 

The BnaGLNl gene expression was monitored at the veg- 
etative stage measuring transcript levels by quantitative real- 
time RT-PCR in samples of taproot, crown, limbs, and veins 
of plants grown under low or high nitrate conditions. Plants 
grown under low or high nitrate conditions grew 13 and 17 
leaves, respectively. F Y /F m and SPAD measurements on all the 
leaf ranks (numbered from the bottom leaf to the top leaf) 
were done to estimate the relative leaf senescence status of 
each leaf. From both SPAD and F Y /F m as senescence markers, 
six leaves were selected from each nitrate condition presenting 
differential senescence levels to perform further experiments 
(Fig. 4). Leaves of rank 3, 5, 6, 7, 9, and 1 1 were harvested on 
plants grown under low nitrate conditions. Leaves 3, 5, 6, 9, 
12, and 15 were harvested on plants grown under high nitrate 
conditions. To simplify the presentation of further results, 
the collected leaf ranks were renamed 1, 2, 3, 4, 5, and 6, 
respectively, with 1 designating the bottom-most and oldest 
collected leaf and 6 the youngest collected leaf. Leaves dis- 
sected as limbs, and primary and secondary veins were used 
to measure BnaGLNl.l gene expression levels in the different 
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Limbs 



2 nd veins 



Fig. 4. Leaf senescence markers on vegetative B. napus plants. F v /F m photosystem II capacity (A, B) and chlorophyll relative content (C, D) were 
monitored at the vegetative stage in all the leaf ranks of four B. napus plants grown under low (white bars) or high (black bars) nitrate conditions. The 
expression of the BnaA.GLN2.a and BnaC.GLN2.a marker genes was quantified (E, F) on selected leaf ranks (*) and confirmed differential senescence 
symptoms. Mean and standard deviation of four plant repeats are shown. 



tissues. In addition to BnaGLNl.l expression, the expression 
of BnaGSLl and BnaGSL2 encoding the chloroplastic GS2 
isoenzymes was also monitored and used as a control for leaf 
senescence as it is known that genes encoding GS2 izoen- 
zymes are down-regulated with leaf ageing in all the plant 
species studied so far (Masclaux-Daubresse et al, 2008). 
BnaGSLl and BnaGSL2 expression levels confirmed the 



differential senescence phenotype of the chosen leaf ranks. 
Leaves 1, 2, and 3 can be considered as senescing leaves, 4 
and 5 as mature leaves, and 6 as a young leaf according to 
Masclaux et al (2000) (Fig. 4E, F). 

Genes that are preferentially expressed under high or low 
nitrate conditions were identified. The results showed that 
regarding the N regime, all the members of the same gene family 
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respond similarly, with a few exceptions from the BnaGLNl.3 
family. It was observed that all the members of the BnaGLNl.l 
and BnaGLNl.4 gene families were significantly induced under 
low nitrate conditions compared with high nitrate. This was 
observed in limbs, secondary veins (Fig. 5A-F), and also in 
primary veins for some BnaGLNl.4 genes (Supplementary Fig. 
S5E-H at JXB online). Induction under low nitrate conditions 
was also clearly observed in the taproot and crown (Fig. 6A, 
B; E-H). In contrast, the two BnaGLN1.2 genes were signifi- 
cantly more expressed under high nitrate conditions in leaf 
limbs and veins but not in the taproot and crown (Fig. 5M, N; 
C, D). Finally no difference was observed in the expression of 
the BnaGLNl.3 and BnaGLN1.5 families in leaf limbs or veins 
between the high nitrate and low nitrate conditions (Fig. 5G- 
L). Surprisingly, all the BnaGLNl.3 genes are significantly 
induced under low nitrate conditions in the taproot but not in 
crown tissue (Fig. 6I-M). Therefore, N-dependent regulation 
might be different in the root and shoot. 

BnaGLNl genes also appeared to be differentially 
expressed depending on leaf ageing and senescence. The two 
BnaGLNl.l genes and the four BnaGLNl.4 genes were signif- 
icantly induced in leaf limbs and veins with ageing and during 
senescence independently of nitrate conditions (Fig. 5A-F; 
Supplementary Fig, S5B, E-H at JXB online). Cumulated 
expression of the BnaA.GLN1.5.alBnaC.GLN1.5.a genes 
(Fig. 5L) also increased with leaf ageing and senescence in 
limbs and veins of plants grown under low and high nitrate 
conditions. In contrast, senescence triggers an opposite effect 
on the mRNA level of the two BnaGLNl. 2 genes, especially 
under low nitrate nutrition (Fig. 5M, N). The effect of senes- 
cence was less evident under high nitrate conditions in limbs 
and veins due to the already high BnaGNL1.2 expression in 
mature leaves. In these leaves, expression profiles are bipha- 
sic, increasing from young to mature leaves then decreasing in 
senescing leaves. Profiles are more complex in the BnaGLNl.3 
family since the expression of BnaC.GLNU.a, BnaA. 
GLNU.b, and BnaC. GLNU.b (Fig. 5H-J) is repressed with 
leaf ageing in limbs and veins, while the expression of BnaA. 
GLN1.3.a (Fig. 5G) and cumulated BnaA. GLN1.3. clBnaC. 
GLN1.3.C (Fig. 5K) is increased with ageing in limbs. 

These results show that within all the BnaGLNl families 
except BnaGLN1.3, all members show similar expression lev- 
els. The four BnaGLNl.4 genes are the most highly expressed 
in all the tissues studied. BnaGLNl.4 gene expression is four 
times higher than that of the BnaGLNl.l genes and 20 times 
higher than that of BnaGLNl. 2. The expression level of all 
the BnaGLNl.3 and BnaGLNl. 5 genes is much lower, except 
that of the cumulated BnaGLNl. 3c genes that reach a similar 
level to BnaGLNl. 2. Table 5 summarizes the N and senes- 
cence effects observed on the BnaGLNl expression levels. 

BnaGLNl genes are differentially expressed at the 
reproductive stage depending on plant organs or 
leaf ageing 

In order to monitor BnaGLNl gene expression at the repro- 
ductive stage, plants were grown in field conditions under low 
or high N regimes. Two leaf ranks (young top leaf and old 



bottom leaf) and the two corresponding stem sections (also 
referred to hereafter as young and old) were collected at flow- 
ering and during grain filling. 

Globally, effects of N limitation on BnaGLNl expression 
were similar to those found at the vegetative stage, except that 
the magnitude of gene repression or induction was lower than 
that observed at the vegetative stage (Supplementary Table S5 
at JXB online). 

Figure 7 reports the effect of senescence on the expression 
of the BnaGLNl genes in leaves and stems of plants grown 
under a sufficient N regime. As a control of leaf senescence 
stages, the BnaGSLl and BnaGSL2 genes are significantly 
more highly expressed in the young tissues than in older tis- 
sues (Fig. 7N, O). There was a sharp decrease in BnaGSL 
gene expression at the flowering stage, while at the seed filling 
stage the magnitude of BnaGSL repression was much lower 
but still significant. 

As observed at the vegetative stage, the BnaGLNl.l genes 
were up regulated with leaf and stem senescence, but this was 
only observed at the flowering stage (Fig. 7A, B). During seed 
filling, expression in leaves and stems was higher than dur- 
ing flowering, showing an effect of plant ageing. However, no 
difference was observed between the young and old leaves, 
suggesting that both types of leaves had become senescent 
between flowering and seed filling. The two BnaGLNl.l 
genes were expressed more highly in leaf blades than in stems 
at both the flowering and seed filling stages. Similarly the 
BnaGLNl. 2 genes were more expressed in leaf blades than in 
stems (Fig. 7C, D). The effect of senescence on BnaGLNl. 2 
genes was opposite to the effect observed on BnaGLNl.l 
genes. BnaGLNl. 2 expression decreased 2- to 3-fold in old 
leaf blades and old stems compared with young leaf blades 
and young stems, respectively. The biphasic profile obtained 
for BnaA.GLN1.2 at the vegetative stage was also observed 
at the flowering stage (data not shown). As observed with the 
BnaGLNl.l genes, the effect of senescence was no more sig- 
nificant at seed filling. 

Among the four BnaGLNl.4 genes, only BnaA.GLNlA.a 
and BnaC.GLNlA.a shared similar expression profiles 
(Fig. 7E-H). They are preferentially expressed in leaf blades 
rather than in stems. In contrast to the vegetative stage, 
BnaA. GLNlA.a and BnaC. GLNlA.a tend to be repressed by 
senescence in leaf blades but induced by senescence in stems. 
This trend was especially significant at the flowering stage. In 
contrast, BnaA.GLNlA.b was induced by senescence in leaf 
blades and stems at the flowering stage but repressed during 
seed filling (Fig. 7G). Finally, BnaC.GLNlA.b expression 
was higher in leaf blades than in stems and was repressed by 
senescence at the flowering stage, similarly to the two BnaA. 
GLNlA.a and BnaC.GLNlA.a genes (Fig. 7H). In contrast 
to the vegetative stage, the members of the BnaGLNl.4 fam- 
ily have developed specificities and are differentially expressed 
at the flowering and seed filling stages. It is likely that they 
have different roles and influences on N metabolism after 
flowering. 

Among the BnaGLNl.3 members, similar profiles 
were observed for BnaA.GLNU.a, BnaA. GLNU.b, and 
BnaC.GLN1.3.b (Fig. 71, K, L). These three genes are 



3940 I Orselefa/. 



30 
20 
10 
0 



c 

O 140 

"55 120 

(/) 100 



BnaA.GLN1.1.a 



I 



Li 



Q. 

X 
<D 

0) 
> 

a: 



to 
'E 

3 

c 
o 

"35 

(/> 
CD 

L- 

Q. 

X 

a> 

o 
> 



80 
60 
40 
20 
0 
50 
40 
30 
20 
10 
0 

5 
4 
3 
2 
1 
0 

0.3 

0.2 



C BnaA.GLN1.4.a 



ILL 



Et T 



ilLiLl 



30 
20 
10 



B 



BnaC.GLN1.1.a 



ilk 



5. FLi. jka. » . 2 SM 



140 I D BnaC.GLN1.4.a 

120 - ~~ 

100 - 

80 - 

60 - 

40 - 



BnaA.GLN1.4.b 



i 



u 



50 
40 
30 
20 
10 
0 



F 

Hi J 



i 



i 



1 



I 



LIMi 



BnaA.GLNI. 3. a 



5 -i T H BnaC.GLN1.3.a 
4 ■ 



3 ■ 
2 - 
1 - 



as 

O 0.1 

0 
20 

(0 15 

"| 10 

§ 5 

3 0 

<» £ 
O 

> f 

JS 4 " 
0) 

a. 2 h 



BnaA.GLNI .3.b 



1 



id 



o 

0.3 
0.2 
0.1 

0 



ikiilikil 



BnaC.GLN1.3.b 



n-.3Ln_.fi. jt 



i 



K 



BnaA.GLNI. 3.c/BnaC.GLN1.3.c 0.02 n 



0.01 - 



BnaA.GLNI. 5.a/BnaC.GLN1.5.a 



M 



BnaA.GLNI. 2.a 



10 



N I hi 1 1 iin 



III. 


1.1a 






■ N 

: iwi 


BnaC.GL 

m 





123456123456 



123456123456 



Limbs 



2 nd veins 



Limbs 



2 nd veins 



Fig. 5. Expression of BnaGL.N1 genes is modified depending on nitrate availability and leaf ageing. The relative expression level of BnaGL.N1 genes 
was monitored in limbs and secondary veins of six leaf ranks harvested on vegetative plants grown under low (white bars) or high (black bars) nitrate 
conditions. Leaf ranks represented as number 1 (bottom and older leaf) to 6 (top and younger leaf) showed differential senescence symptoms. Mean and 
standard deviation of four plant repeats are shown. 



Exploring glutamine synthetase gene expression in Brassica napus L | 3941 



2.5 

2 

1.5 
1 

0.5 
0 



BnaA.GLNl.l.a 



ill 



BnaC.GLNl.l.a 



ill 



2.5 
2 

1.5 
1 

0.5 



BnaA.GLN1.2.a 



BnaC.GLN1.2.o 



2.5 



1.5 



0.5 



El 




BnaA.GLNlA.a 



BnaC.GLNlA.a 



5 i 

4 
3 
2 
1 
0 



i 



BnaA.GLNlA.b 

0.8 -i G 0.2 



t 



0.6 



0.4 



0.2 



1 0.1 
-i 1 ^ i 



BnaC.GLNlA.b 



BnaA.GLN1.3.a 



BnaC.GLN1.3.a 



2.5 
2 

1.5 
1 

0.5 
0 

40 
30 
20 
10 



1 



M 



BnaA.GLNl.3x/ 
BnaC.GLN1.3.c 



k 



Taproot 



Crown 



20 



15 



10 



BnaA.GLN1.3.b BnaC.GLN1.3.b 
K 0.3 n L 



uu 



It 



0.006 
0.005 
0.004 
0.003 
0.002 
0.001 
0 



BnaA.GLN1.5.a/ 
BnaC.GLN1.5.a 5 



ill 



Taproot 




Crown 



Taproot 



Crown 



Taproot 



Crown 



Fig. 6. Expression of BnaGLNI genes is modified depending on nitrate availability in the taproot and crown of vegetative B. napus plants. The relative 
expression level of BnaGLNI genes was monitored at the vegetative stage in the taproot and crown of four plants grown under low (white bars) or high 
(black bars) nitrate conditions. Mean and standard deviation of four plant repeats are shown. 



down-regulated in old leaves and stems compared with 
young leaves at the flowering stage. However, their expression 
increased sharply in old limbs at the seed filling stage. The 
other BnaC. GLNlJ.a and BnaA. GLNlJ.cl BnaC. GLN1.3.C 
expression profiles did not show any modification associated 
with leaf or stem senescence (Fig. 7 J, M). All the BnaGLNI. 3 
genes appeared to be more highly expressed in stems than in 
leaves, especially at the flowering stage. 

It was not possible to measure BnaA.GLN1.5.clBnaC. 
GLN1.5.C gene expression, possibly due to the very low 
expression level in vegetative tissues that cannot be accurately 
measured in field-grown plants. 

Results obtained at the flowering and seed filling stages 
confirm results from the vegetative stage. BnaGLNI genes 



are generally similarly regulated according to their orthology 
group, although exceptions were observed particularly at the 
seed filling stage, such as with BnaA.GLNlA.b and BnaC. 
GLN1.3.a(¥ig. 7G, J). 

Discussion 

Glutamine synthetase is a key enzyme of N metabolism 
involved in ammonium assimilation and remobilization. 
Recent studies highlight the important role of GS1 cytosolic 
isoenzymes for N management linked to yield establish- 
ment and seed filling in monocotyledonous crops (Tabuchi 
et al., 2005; Martin et al., 2006; Bernard and Habash, 2009; 
Swarbreck et al, 2011). The GSl-coding genes are therefore 
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Fig. 7. BnaGLNI genes are differentially expressed depending on flowering or seed filling stages and leaf age. The relative expression level oi BnaGLNI 
genes was monitored in young and old leaves of plants grown in the field and supplemented with nitrogen. Leaf limbs and stems were collected at 
flowering (white bars) and seed filling (black bars) stages. Mean and standard deviation of four plant repeats are shown. * indicates significant difference 
(Student's f-test, P<0.05) between flowering and seed filling stages. 



good candidates for improving yield traits and grain quality 
(Masclaux-Daubresse et al, 2008). The complexity of study- 
ing glutamine synthetases arises from the fact that two isoen- 
zymes exist, one in the chloroplast and the other in the cytosol, 
and that several isoforms exist for the cytosolic enzyme. The 
numerous isoforms are encoded by a multigenic family, and 
the five GLN1 genes in A. thaliana are likely to present dif- 
ferent roles depending on plant organs and nitrate availabil- 
ity in the soil (Lothier et al, 2011). Similarly the five maize 
GLN1 genes do not participate equally in N management at 
the whole-plant level (Martin et al, 2006; Hirel et al, 2007). 
The aim of this study was to identify the whole BnaGLNI 
gene family and to characterize the expression of the different 
genes depending on nitrate availability as well as depending 
on ageing and leaf and stem tissue senescence. 

Using the sequences obtained from EST libraries and 
genome sequencing, a total of 16 genes belonging to the 
BnaGLNI family, eight genes from each of the A and C 
genomes, were found. In accordance with the history of 
B. napus genome formation (Nagaharu, 1935), it was found 
that each BnaGLNI gene is closely related to a BraA. GLN1 
or BolC. GLN1 gene depending on its A or C genome loca- 
tion. Therefore, it can be stated with confidence that all 
the GLN1.1 genes of B. napus have been described in this 
report. Sequence analyses also showed that B. napus coding 



sequences are highly conserved between the A and C genomes 
and also between B. napus genotypes. The level of sequence 
divergence observed in the BnaGLNI family is between 
0.9% and 2.9% SNPs (single nucleotide polymorphisms) in 
CDS, which is less than the preliminary observation show- 
ing sequence divergences of -3-5% SNPs in CDS from SLR1 
(Inaba and Nishio, 2002). 

It is well known that B. napus shows a high degree of collin- 
earity to its diploid progenitors B. rapa and B. oleracea (Rana 
et al, 2004). Many studies have investigated the segmental 
structure of the Brassica genomes and led to the identification 
and genetic mapping of syntenics blocks between A. thaliana 
and the Brassica genomes (Parkin et al, 2005; Schranz et al, 
2006; J. Wang et al , 201 1). The number of potential BnaGLNI 
genes and their localization on linkage groups can then be 
predicted depending on the number of times the blocks are 
replicated and on their localization on each B. napus link- 
age group. According to the whole-genome triplication event 
that occurred in Brassicaceae genome species after divergence 
from Arabidopsis (Lysak et al, 2005; Parkin et al, 2005), and 
the recent hybridization between B. rapa and B. oleracea lead- 
ing to the appearance of B. napus, each AtGLNl gene could 
have been found in triplicate in each A and C genome from 
B. napus to form three pairs of homeologous genes. From 
EST and genome sequence analyses, it is revealed here that 



3944 I Orselefa/. 



only one homeology group exists for AtGLNl.l, AtGLNl. 2, 
and AtGLN1.5 and two groups for AtGLNl.4. AtGLN1.3 is 
the only GLN1 gene for which the six BnaGLNl orthologous 
genes were retained in the B. napus genome. This illustrates the 
massive gene loss that occurred in the Brassica lineage after 
the whole-genome triplication event. Indeed, using A. thaliana 
as an outgroup, Town et al. (2006) found that 35% of genes 
inferred to be present when genome triplication occurred in 
the Brassica lineage have been lost in B. oleracea. Similarly, 
whole-genome analysis of B. rapa revealed a high rate of gene 
loss, from 30% to 64% depending on the degree of fractiona- 
tion of the region considered (X. Wang et al , 201 1). BnaGLNl 
families are a good example of this. 

With the exception of BnaGLNl. 4.b genes, the present study 
points out that pairs of BnaGLNl homeologous genes share 
very similar transcription profiles. Furthermore, within one 
orthology group, when several groups of homeologues were 
retained, paralogous genes conserved similar expression pro- 
files (BnaGLNl. 3 and BnaGLNl. 4). This suggests that coding 
but also regulatory sequences were essentially conserved after 
genome merging of B. napus progenitors, but also after the 
whole-genome duplication (WGD) and diploidization events 
that occurred in the Brassica lineage after the divergence from 
the Arabidopsis genus. WGD is generally thought to provide 
raw material for gene neo- and subfunctionalization, extending 
resilience to deleterious mutations, increasing the net speciation 
rate and species richness (Soltis et al, 2009), as well as provid- 
ing the adaptive advantage for colonizing harsh and unstable 
environments (Franzke et al, 2011). On the other hand, main- 
tenance of redundancy can confer robustness against mutations 
(De Smet and Van de Peer, 2012) and/or a selective advantage 
in increasing the abundance of encoded proteins (Bekaert et al. , 
201 1). As GS1 is an essential enzyme of primary N metabolism, 
linked to central carbon metabolism via the GS/GOGAT cycle 
that might also play a role in the adaptation of plant to nutrient 
deficiency and pathogen attack (Brauc et al, 2011; Seifi et al, 
2013), maintenance of multicopy of GLN1 could confer robust- 
ness against mutations. Nevertheless, GLN1 expression profiles 
have not been exhaustively investigated and there might be par- 
ticular environmental or developmental conditions allowing 
the differentiation of expression profiles between homeologous 
and/or paralogous genes. Partially overlapping profiles could 
provide robustness against mutations, but also adaptive advan- 
tages for colonizing harsh and unstable environments. 

Allopolyploidization involves the merger of two different, 
and often divergent, genomes whose reconciliation in a com- 
mon nucleus often leads to myriad changes, including unequal 
expression of the two merging genomes. Biased expression 
among homeologues has been found in cotton and wheat 
(Pumphrey et al, 2009; Rapp et al, 2009). Previous studies 
have suggested bias toward the B. rapa A genome in the tran- 
scriptional expression of rRNA genes (Chen and Pikaard, 
1997). In the present study it is shown that most of the 
BnaGLNl homeologous pairs display similar expression levels 
in the various tissues studied. Differences in mRNA contents 
observed between homeologous pairs of A and C genome ori- 
gins were generally very small. No systematic bias towards the 
B. oleracea C parental genome or B. rapa A genome can be 



identified in this study. Bias towards the B. oleracea C parental 
genome was identified for BnaGSL (at flowering and seed fill- 
ing stages), BnaGLNl. 3. a (especially at the vegetative stage), 
and BnaGLNl. 4. a homeologues, and in favour of the B. rapa 
parental genome for BnaGLNl. 1. a and BnaGLNl. 2. a. These 
results are in agreement with the recent finding that nearly 7% 
of the potentially identified homeologous genes expressed in 
a leaf extract are displaying a differential expression level in 
favour of the A or C parental genome for 1/3 and 2/3 of the 
pairs, respectively, and that genes involved in metabolic pro- 
cessees tend to be over-represented (Higgins et al, 2012). 

Beside their intrinsic expression levels, it was found that 
the BnaGLNl genes are similarly regulated depending on 
their orthology group and that they are differentially regu- 
lated between groups. Overall it was found that specificities 
of expression are conserved between BnaGLNl genes and 
their respective AtGLNl orthologues, raising the hypothesis 
of conserved physiological functions. 

In Arabidopsis, several studies have shown that AtGLNl.l 
is highly induced in leaves during senescence (Guo et al, 
2004) and is up-regulated when exogenous N sources are 
limiting (Lothier et al, 2011). Up-regulation of AtGLNl.l 
under nitrate starvation is in good agreement with the kinetic 
properties described by Ishiyama et al. (2004) that suggested 
that the high affinity of AtGLNl.l for ammonium is cor- 
related with a role for the enzyme under low N conditions. 
The induction of the expression of the two BnaGLNl. 1 genes 
in older stems and leaves is conserved, as already shown by 
Buchanan- Wollaston and Ainsworth (1997) and Ochs et al. 
(1999). It was also found that the two BnaGLNl. 1 genes 
are overexpressed under low nitrate conditions especially at 
the vegetative stage. As the amino acid residues known to 
be involved in the high affinity of AtGLNl towards ammo- 
nium are not conserved in any of the BnaGLNl proteins, 
no information about the potential kinetic properties can be 
extrapolated from the protein sequence. Regarding expres- 
sion profiles, it is at least suspected that BnaGLNl. 1 and 
AtGLNl.l proteins might have similar roles 

In a previous study, it was found that AtGLNl. 2 is slightly 
induced by leaf ageing and that the expression profile is 
biphasic, with an increase from young to mature leaves and 
then a decrease in strongly senescing leaves (Diaz et al, 2008; 
Lothier et al. , 201 1). In addition, it was found that AtGLNl. 2 
is mainly expressed in roots and leaves under a high N regime 
(Guo et al, 2004; Ishiyama et al, 2004; Lothier et al, 2011). 
A detailed functional analysis led to the conclusion that 
AtGLNl. 2 was involved in primary ammonium assimila- 
tion under high N regimes (Ishiyama et al, 2004; Lothier 
et al, 2011). BnaGLNl. 2 genes (also named BnGSR2.1 and 
BnGSR2.2) are also more highly expressed in roots than 
in shoots (Ochs et al, 1999). It is shown here that the two 
BnaGLNl. 2 genes are more expressed in young leaves than in 
old leaves and are overexpressed under high N regimes, sug- 
gesting a similar role to AtGLNl. 2 in primary ammonium 
assimilation. Interestingly, a BraA.GLN1.2 gene, named 
BcGSl (Sun et al, 2010), was found also to be expressed in 
root and induced under high N regimes, suggesting the con- 
servation of the regulation in the Brassiceae tribe. 
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Similarly to AtGLNl. 3, the BnaGLNl.3 genes are not 
induced in older leaves at the vegetative stage and do not 
respond to differential N regimes (Guo et al , 2004; Lothier 
et al, 2011). The hypothesis about the physiological roles of 
AtGLN1.3 in N export via phloem tissues in roots is sup- 
ported by the high capacity of the enzyme for glutamine 
synthesis and by the location of AtGLN1.3 expression in the 
root vasculature (Ishiyama et al, 2004). It has to be noted 
that BnaGLNl.3 genes are the only BnaGLNl genes prefer- 
entially expressed in stem tissues compared with leaf blades. 
As stems are richer in vascular tissues than leaf blades, the 
higher expression of BnaGLNl.3 in stems might be related to 
a potential vascular localization that remains to be explored. 

AtGLNl.4 is one of the markers used for leaf senes- 
cence (Guo et al, 2004; Wingler et al, 2009). AtGLNl.4 is 
induced by N limitation or starvation in both the root and 
shoot (Ishiyama et al , 2004; Lothier et al , 201 1). AtGLNl .4 
protein exhibits high affinity towards ammonium and is 
expressed in the pericycle cells of roots (Ishiyama et al, 
2004). Evidence for BnaGLNl. 4a root expression was found 
in EST libraries. However, induction of gene expression in 
leaves in response to ageing and low N regime is not well con- 
served among the four BnaGLNl. 4 orthologues in regards 
to developmental stages. The BnaGLNl. 4 genes are signifi- 
cantly overexpressed in senescing leaves and under low nitrate 
conditions at the vegetative stage. However, this trend is not 
conserved at the flowering and seed filling stages. BnaA. 
GLNlA.a, BnaC.GLNlA.a, and BnaC.GLNlA.b are clearly 
and significantly less highly expressed in old than in young 
leaves at these stages. Furthermore, as the residues confer- 
ring the high affinity property are partially conserved in 
BnaGLNl. 4 proteins, it might be suspected that BnaGLNl. 4 
genes have a role when N resources are low. Surprisingly, in 
contrast to Arabidopsis in which AtGLNl.4 expression is one 
of the lowest, BnaGLNl. 4 expression levels are the highest 
found among all the BnaGLNl genes. 

In Arabidopsis, AtGLNl. 5 expression is known as the low- 
est of the AtGLNl gene family. Expression was mainly found 
in seeds (Lothier et al, 2011), and very little is known about 
AtGLNl. 5. BnaGLNl. 5 gene expression is also very low, 
and it was not possible to measure it in leaf and stem tissues 
at flowering and vegetative stages. Similarly to AtGLNl. 5, 
BnaGLNl. 5 ESTs were found in reproductive tissues and 
seeds, thus suggesting specific roles during seed maturation. 

Conclusion 

A total of 16 B. napus GLN1 genes were identified, among 
which 12 have never been described. The total number of 
BnaGLNl genes, their phylogenetic relationships, and genetic 
location are in agreement with the evolutionary history of 
Brassica species. Some specificities of expression seemed to be 
conserved among the Brassiceae tribe and especially between 
A. thaliana and B. napus. Regulations arising from plants 
interactions with their environment (such as N resources), 
final architecture, and therefore sink-source relations in 
planta, seem to be globally conserved when compared with 
data available from the Arabidopsis model. Considering the 



architectural, size, and lifespan differences between A. thali- 
ana and B. napus, it is not surprising to find some differences 
in gene expression profiles. Also, due to the higher number 
of GLN1 genes conserved in the B. napus genome, it seems 
correct to find some specificities in the expression of each 
BnaGLNl in contrast to genes involved in flavonoid bio- 
synthesis that display highly conserved expression profiles 
between A. thaliana and B. napus during seed development 
and are highly dependent on tissue differentiation (Auger 
et al, 2009). A more detailed localization of BnaGLNl gene 
expression would refine the hypothesis concerning their physi- 
ological role. Indeed, the present expression study relied on 
leaf blade and stem samples consisting of different tissues 
with contrasting physiological roles, in particular parenchyma 
and vascular tissues. Such a localization study could be advan- 
tageously performed in Arabidopsis and B. napus to provide a 
new basis for comparison of the evolution of this gene family. 
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Supplementary data are available at JXB online 

Figure SI. BnaGLNl gene localization on the A or C 

genome using a panel of various Brassica genotypes. 

Figure S2. BnaGLNl gene localization on LGs using 

mono- and polysomic additional lines. 

Figure S3. Positions of BnaGLNl genes on the B. napus 

genetic map. 

Figure S4. Alignments of the genomic sequence, deduced 
mRNA sequence, and contig of EST sequences are reported 
for each BnaGLNl gene. 

Figure S5. Expression of BnaGLNl genes is modified 
depending on nitrate availability and ageing in primary veins 
of vegetative B. napus plants. 

Table S 1 . Clones from Genoplante and ADIS-MPIZ oilseed 
rape cDNA libraries used for BnaGLNl mRNA sequencing. 

Table S2. Primers used for cloning and genetic mapping of 
the BnaGLNl gene. 

Table S3. Specific primers used for sequencing 
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Table S4. qPCR primers used for BnaGLNl gene expres- 
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Table S5. Induction of BnaGLNl gene expression under 
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belonging to each B. napus, B. oleracea. and B. rapa GLN1 
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ies according to their BnaGLNl contigs. 
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